Effect of Level Statistics on Superconductivity in Ultrasmall 

Metallic Grains 



Robert A. Smith 1,2 and Vinay Ambegaokar 2 
1 School of Physics and Space Research, University of Birmingham, Edgbaston, 
Birmingham B15 2TT, United Kingdom 
2 Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca NY14853, USA 

(August 26, 1996) 

Abstract 

We examine the destruction of superconducting pairing in metallic grains 

as their size is decreased for both even and odd numbers of electrons. This 

occurs when the average level spacing d is of the same order as the BCS order 

parameter A. The energy levels of these grains are randomly distributed 

according to random matrix theory, and we must work statistically. We find 

that the average value of the critical level spacing is larger than that for the 

model of equally spaced levels for both parities, and derive numerically the 

probability densities P 0je (d) that a grain of mean level spacing d shows pairing. 
PACS numbers: 74.80. Bj, 71.23.-k, 73.20.Dx 
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A recent experiment by Black, Ralph and Tinkham U involving the observation of a 
superconducting gap in ultrasmall Al grains (of size between 5 and 13nm) has led to recon- 
sideration of an old but fundamental theoretical question - how small can a superconductor 
be? It is also of interest that in a previous experiment on a smaller Al grain, the same 
group saw no sign of a gap (although, as they noted, there are experimental difficulties in 
observing an energy gap of similar magnitude to the average level spacing), von Delft et 
al. H have considered a simple mean field model for such a grain which uses the standard 
BCS pairing interaction and assumes equal level spacing for analytical simplicity. Even-odd 
parity effects which can be seen in samples at least 10 4 times larger ||[7j , and are 

of paramount importance here, are included in their model. They find that the supercon- 
ducting gap at zero temperature should cease at critical level-spacing d° = 0.89A(0) in odd 
grains, and d e c = 4d° in even grains, where A(0) is the zero-temperature bulk gap. They also 
note that although the mean-field result is subject to several types of correction, it does give 
a criterion for when pair correlations will cease to exist. It is therefore surprising that in 
the data of BRT from sample 4, an odd grain, that a gap is still seen although the sample is 
very close to the odd critical level-spacing. We also note that BRT's data shows no variation 
of the gap with level-spacing d, suggesting that their samples are still on the flat part of the 
A(d) curve. 

In this paper we consider the effect on the mean-field theory of relaxing the condition of 
equal level-spacing. It is by now well-known that the level-spacing in small metallic grains 
is the Wigner-Dyson (WD) distribution || obtained from random matrix theory (RMT) || . 



This was first conjectured by Gor'kov and Eliashberg [|10j, and later proved by Efetov ||11|| . 
The reason for considering this effect is that most of the other corrections to mean-field 
theory seem to lead to a reduction in d c \ on the other hand, level statistics effects lead to 
larger values of (d c ), as we shall see. 

The first thing we shall do is to reproduce the results of von Delft et al. for d° c and 
d e c . We do this to demonstrate how the positioning of the energy levels enters into the 
calculation, and how this leads to the factor of 4 in the result d e c = 4d°. Our starting point 



is the mean-field self-consistency equation 



\ = d £ ^E.i 1 - 2 ^ (!) 
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where Et = yj (q — fi) 2 + A 2 , with e« the i-th energy level, \i the chemical potential, uj c the 
Debye energy, d the level spacing, and A the BCS interaction. The occupation factor fa 
differs for even or odd parity ensembles 

where = ±(e^ Ei ± and Z± = ± e _/3£i ). We will work at zero temperature, 
so that fi = 1/2 if the chemical potential lies on a level, and zero otherwise. In the case 
of equal level spacing the chemical potential lies half-way between the last filled and first 
empty levels in the even case, and on the half-filled level in the odd shown in Fig. 

(1). For the case of the critical level spacing, the solution has A(T = 0) = 0, so that one 
has 

A = g iTTJl ; A = g f (3) 
These can be rewritten in terms of the digamma function to yield 



j = 4>(u c /dt) - ^(1/2) » log (u c /d e c ) - ^(1/2) 
i = ^(w c /dg) - V(l) ~ log (oj c /d°) - tP(1). 



(4) 



Finally, since we know that = —7, if) (1/2) = —7 — 2 In 2, where 7 is the Euler- 

Mascheroni constant, it follows that 

d° = Vc^e" 1 ^ = ie 7 A(0) « 0.89A(0) ; d e c = 4d° c . (5) 

We see that the factor of 4 between and d° comes from the fact that ip(l) —■0(1/2) = 2 In 2, 
and thus ultimately from the positioning of the chemical potential relative to the energy 
levels. Furthermore if we write this out as a series for 2 In 2, 

2 In 2 = -0(1) -"0(1/2) 

(6) 
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we see that only a few terms are needed before we are close to the final answer. In other 
words, the value of this particular factor is determined by a few energy levels near to the 
Fermi surface. 

Let us now see how we can extend the above approach to the case of non-equally spaced 
energy levels. Consider first the odd case, so that the chemical potential lies on the half-filled 
level eo- The gap equation can then be written as 



We can then take the average of this equation over the disorder ensemble, so that all the 
statistical information about the level spacing occurs in the average over the sum of delta 
functions. This can be related to the two-level correlation function (TLCF) of the system 
as follows. The TLCF is defined by 

R(e-e') = d 2 (j:6(e-e l )5(e'-e J )}. (8) 

Since this is a function only of the energy difference, e — e', we may set e' = to obtain the 
result 



= d a E«(e-e i + c J )%)> (9) 
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The averaged odd-gap equation can finally be written as 

1 fTTUJc/d d X 



, =R(x) (10) 

where x = irui/d. Since the system has time-reversal invariance, the TLCF given by RMT 
is that for the orthogonal ensemble ||, 



and the average critical level spacing (d°) is then the solution of 

1 fTTWc/rf 



A jo x 
This integral can then be performed analytically to give 



(12) 



(d°) = 7re 7+ " 2/16 - 7/4 A(0) = 1.80A(0) (13) 

and we see that the average odd critical level spacing is a factor 2 larger than in equal level 
spacing model. 

Next let us consider the even case. The chemical potential is now halfway between the 
last filled level, eo, and the first filled level, e\. It follows that the gap equation can now be 
written in the form 

\ = dJ2^= 1 (14) 
A Y 2y/(ei - ^) 2 + A 2 

where the sum over i includes both eo and e\. For e^ not equal to eo or ei, we can rewrite 

Q Y~ = (<* - d) + ■ (15) 

We know from the odd case that the distribution of the e« — e\ is described by the TLCF, 
R(x). The distribution of the e\ — eo is given by the nearest level spacing distribution, 
P(y). There is no analytic expression for P(y), but it is well approximated by the "Wigner 
surmise" ||, 

Ply) = ^ y2/in ( 16 ) 

where y = 7r(ei — eo)/d. Let us now assume that the distributions of — e\ and ei — eo may 
be treated independently. The actual distribution function we need is a three-level function 
for eo, ei, and e^, but such a function is not discussed in the RMT literature. The equation 
for the critical level spacing, d% is then 

1 f°° , r./ x 2?r r™*/ d c , f°° , R(x)P(y) 

-=/ dyP(y)—+ dx dy K > ) yj 17 
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where the first term comes from the levels eo and ei, which have to be treated separately, and 
the second term comes from all other levels. Note that if we were to replace the denominator 
x + y/2 by x in the second term we should recover the odd integral. It follows that we should 
evaluate the difference between the second term and the odd integral, from which we obtain 
the result 

(^>=exp(7r-2//7r)(d2> (18) 

where / is the integral 

roo „ roo I 

/=/ dtt 2 e- t2 /* - x 
Jo Jo x(x + t) 

r • 2 , ■ xl ( 19 ) 
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This integral cannot be performed analytically, and has the numerical value / ~ 1.7343. We 
can therefore summarize the results for the mean critical spacings in terms of the bulk BCS 
gap A(0) or the critical spacing for equidistant levels in odd grains d° by 

(d°) = 1.80A(0) « 2.0d° , (d e c ) = 7.67(d°) w 15.5d°. (20) 

We see that the consideration of level statistics not only makes both the odd and even 
critical level spacings larger, it also increases the ratio between them. The reason for this 
is that both the individual gap equations, and the difference between the gap equations, 
involve the inverse of energy level spacings. The fluctuations to smaller level spacings are 
thus weighted more than those to larger than average level spacings i.e. (1/5E) > 1/{SE). 

In the analytic discussions above we have evaluated the mean value of the critical level 
spacing. We note that the mean is only one statistical measure of a probability distribution, 
and may not actually be the one we want. We would therefore like to look at the probability 
distributions P ^ e (d) of there being a superconducting gap in odd/even grains with average 
level spacing d. We might imagine an experiment in which many grains of the same nominal 
size are produced and examined for the presence of a superconducting gap; the experimental 
results would then yield P efi (d). To obtain P etQ (d) we proceed numerically, obtaining sets of 
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energy levels {e{\ by diagonalizing N x N random matrices. Since the eigenvalues produced 
have a semicircular density of states ||, 

1 



p( e ) = -V2iV - e 2 9{2N - e 2 ) (21) 

7T 

where 9{x) is the Heaviside function, we use the rescaling 

e -> — [2JV sin^ 1 (e/v 7 ^) + ev^iV - e 2 l (22) 

27T L - 1 

to obtain eigenvalues with average spacing unity |12| . From the gap equations we see that 
the criterion for a grain to have a non-zero superconducting gap is (all energies now in units 
of d) 

\ K E 2^ (23) 
where /i = e in the odd case, and /x = (e + ei)/2 in the even case. As in the analytical 
calculations, d enters only in the upper cut-off u c /d. We choose the value A = 0.193 
corresponding to u) c /d° = 100 in the equal level spacing case, and d is measured in units 
of this d°. For ci-points from 0.5 to 1.0 we use 1000 realizations of 500 x 500 matrices; for 
d-points from 1.0 to 4.0 we use 1000 realizations of 300 x 300 matrices; and for <i-points 
from 4.0 to 20.0 we use 10000 realizations of 100 x 100 matrices. We note that for larger d, 
where we need more realizations to get good statistics, we are fortunate in that we require 
smaller matrices. The results are shown in Fig. (2). 

We see that for both odd and even cases there is significant chance of superconductivity 
persisting beyond the critical level spacings deduced from the equal level spacing model. 
Both curves also show long tails which we believe are due to the non-zero probability of 
finding two levels very close together. These long tails make it hard to estimate the mean 
value of the critical spacing from the numerical data — though it is worth noting that Fig. 
(2) is quite consistent with the analytic results in Eq. (20). The long tails also imply that 
the mean is perhaps not the best measure for a typical critical spacing. If instead we the 
value of d where P(d) = 0.5, we get about l.7d° for the odd case, and 6.5d° for the even 
case. So, using this measure we recover a factor of order 4 between the two critical spacings. 



Let us now discuss the implications of the above calculation for experiment. First let us 
ask the question of what the BRT experiment actually measures for the cases of even and 
odd grains. In both these cases, since an electron tunnels onto and then off the island, the 
result involves some sort of comparison between odd and even states. So do we see the odd 
gap, even gap, or some mixture thereof? We will always work in the zero temperature limit, 
which is effectively where the experiment is performed. Consider the case of the even grain. 
The lowest state an electron can tunnel into is the first unoccupied level, so the energy cost 
should be the energy difference between ground states of the system with 2N and 2N + 1 
electrons. From the T = limit of Eq. (4.9) of Ref. ||, this is given by 

E(2N + 1) — E(2N) = fi + A. (24) 

To see which A is involved, note that A arose from formulas (2.19) of Ref. |4| which give 



fl„ — £l P = — In 



(3 \Z + -Z_J 
l ln fl + (l-2N eff e-^ 



(25) 



(3 \\-{l-2N eff e-^))- 

At zero temperature we see that it is the odd gap, A G , that is measured. 

For the case of an odd grain, the lowest state for the electron to go into is the singly 
occupied state, and we measure the energy difference 

E{2N + 2) - E{2N + 1) = \i - A . (26) 

It follows that the lower branch of the odd grain curve should be the mirror image of the 
even grain curve, so that A D appears in experiments on both even and odd grains. We note 
that the first excited state of an odd grain is obtained by putting the extra electron into 
the first unoccupied state, thus giving a state with two unpaired electrons. This should 
have an energy roughly 2A above the ground state, but the evaluation is complicated by 



the quasi-particles reducing the phase space for pairing correlations. [13 
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FIG. 1. The positioning of the chemical potential relative to the electronic energy levels in a 
superconducting grain with (a) an even number of electrons; (b) an odd number of electrons. In 
the even case the chemical potential lies halfway between the last filled and first empty level; in 
the odd case it lies on the half-filled level. Although illustrated for the equal level spacing case, 
the same occurs for randomly spaced levels. 
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FIG. 2. The probability densities P e ,o(d) of a metallic grain with mean level spacing d, and 
an even (+ symbols) or odd(x symbols) number of electrons, having a non-zero superconducting 
energy gap. d is measured relative to dP c = 0.89A(0), the critical level spacing for the equidistant 
model with an odd number of electrons 
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